Magnetic Bianchi I Universe in Loop Quantum Cosmology 
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We examine the dynamical consequences of homogeneous cosmological magnetic fields in the 
framework of loop quantum cosmology. We show that a big-bounce occurs in a collapsing magnetized 
Bianchi I universe, thus extending the known cases of singularity-avoidance. Previous work has 
shown that perfect fluid Bianchi I universes in loop quantum cosmology avoid the singularity via 
a bounce. The fluid has zero anisotropic stress, and the shear anisotropy in this case is conserved 
through the bounce. By contrast, the magnetic field has nonzero anisotropic stress, and shear 
anisotropy is not conserved through the bounce. After the bounce, the universe enters a classical 
phase. The addition of a dust fluid does not change these results qualitatively. 



I. INTRODUCTION 



Loop quantum cosmology (LQC) is a theory of quantum cosmology based on the more general theory of loop 
quantum gravity [l[ (for reviews, see [1, H, 13 ) . One of the most important predictions of LQC is the avoidance of the 
big-bang singularity, which is replaced by a bouncing universe @, la 0> Hi @] f° r isotropic models sourced by a massless 
scalar field. This result has been derived rigorously at the level of the quantum theory, but has also been understood 
at the level of approximate effective classical equations that capture the main features of the quantum dynamics. 

Extending the rigorous quantum dynamics to the anisotropic Bianchi I model is challenging. Because of this, the 
dynamics of Bianchi I in LQC have mostly been studied by extrapolating the approximate effective semi-classical 
equations that proved successful in the isotropic case. For matter with zero anisotropic stress, the effective equations 
~~ • predict a bounce that avoids the classical singularity [Io| . thus extending the results of the isotropic case. Furthermore, 
it was shown that the shear anisotropy does not blow up in the collapsing phase, but remain finite through the bounce. 
Several ambiguities in the quantum construction of the Hamiltonian constraint were considered [To| . one of which 
has since been favored by gauge considerations [Til ] and by a more thorough construction of the quantum theory [l2| . 
With this choice of quantization scheme, it was shown that the anisotropic shear is in fact conserved across the bounce 
when the matter has zero anisotropic stress [Toj . 

In this paper, we couple a homogeneous magnetic field to a Bianchi I universe and consider the effective semi- 
qq ' classical modifications to the equations of motion. We show that the singularity is still avoided via a bounce, during 
which anisotropies remain finite. However, the anisotropic stress in the magnetic field leads to a non-conservation of 
r* ■ shear anisotropy through the bounce, in contrast to the case where of matter has zero anisotropic stress. 

II. CLASSICAL EQUATIONS 

The inclusion of a cosmological magnetic field breaks isotropy, so we consider an anisotropic Bianchi cosmology, 
the simplest being the Bianchi I model: 

ds 2 = -dt 2 + a 2 1 (t)dx 2 + al{t)dy 2 + al(t)dz 2 , H, := — . (1) 

(Xn 



Maxwell's equations are 



V[»F va] =0, V„i^ = J", (2) 
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where J M is the four-current. The Faraday tensor defines electric and magnetic fields relative to observers with 
four- velocity m m [l3[ 



— Fpi/U , B ^ — ~£}ivaF : (3) 



where e^ va is the alternating tensor in the observer's rest-space. 

We assume high conductivity in the early universe, so that the electric field is effectively zero, and Maxwell's 
equations reduce to 



K"B V = B', (4) 
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bTVpBv = 0, e llva V u B a = h llv J v , (5) 
where /i M „ = g^ v + u^u v projects into the rest-space, O is the volume expansion and a^ v is the shear: 

9 = Hi + H 2 + H 3 = 3- , a 3 := ai a 2 a 3 , of = Sj ( 2* - - ] . (6) 
a \a,j a J 

For a homogeneous magnetic field in Bianchi I, the divergence constraint is automatically satisfied, and the curl 
constraint shows that there is no 3-current. To solve the induction equation 0$, we assume without loss of generality 
that the magnetic field is aligned along the ^-direction: B^ = Bi(t)<5 A1 1 . The solution is 



^=4. ( 7 ) 

{a 2 a 3 y " 6 



where j3 is constant, in agreement with [14j |. 

The electromagnetic energy-momentum tensor is given by 
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where the magnetic energy density and anisotropic stress are 



PBU^Uv + -pBh^ + rf„ , (8) 



Pb = , tt* = ^B a B a h^ u - B^B V . (9) 

For the case with only a magnetic field, the Einstein field equations G^ u — SttGT^ lead to [TJ] : 

Hx+Hj = ^, 7 = 2,3, (10) 

Hfa 6 = l2l3 ~iTrG(3 2 aj, (11) 

where 7/ are constants. It follows that 

7273 > , cti < ai m = ^Qp 2 > ( 12 ) 

and thus any expansion in the magnetic field direction will eventually come to rest at the maximum scale factor ai m 
and turn around into a contracting phase. The solutions a,(f) can be given analytically [bj ]: 



a oc (l + .f)f ±a -\ f(ai):=^ 

ai 



1-4/1- 
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72 +73 n „s 

a := — = , (13) 



V7273 



ai cx (l + / 2 )/±7x/v^W, (14) 

where the ± refers to the expanding (collapsing) branches before (after) a = ai m is reached. 

We can now analyze the singularity behavior of the general solutions. From Eq. (| 1 3(1 . the volume goes to zero or 
infinity when a± goes to zero, depending on the value of a. There are two separate cases. The first is the axisymmetric 
case, a 2 — a 3 , 72 = 73, a = —2. From Eqs. (fT3"|) and (|14[) . the axisymmetric singularity is characterized by 

»2 = »3 — > const , ai,a— >0. (15) 
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Directions orthogonal to the magnetic field freeze as the singularity is approached, while the magnetic field direction 
contracts to zero, along with the total volume. This singularity is present in all solutions for the axisymmetric case. 
The overall evolution is characterized by expansion in the direction of the magnetic field until a\ m is reached. After 
that, the a\ direction contracts and reaches the singularity in finite proper time. For this type of evolution, past 
infinity is characterized by a, a 2 , 03 — > 00 while a\ — > 0. In addition, the time reversed scenario is possible. 

The non-axisymmctric case is slightly more complicated, but again all trajectories are singular. In this case 72 7^ 73, 
and the singularity is characterized by 

a 2 ^0, a 3 — > 00 if I72I > |7s| , a\, a— >0, (16) 

and vice- versa if I72I < I73I. Thus the singularities are again given by a\ going to zero, with one of the orthogonal 
directions contracting to zero, while the other expands to infinity in such a way that a still goes to zero. 



III. EFFECTIVE LOOP QUANTUM EQUATIONS 

The loop quantum formulation is based on a Hamiltonian framework where the gravitational degrees of freedom in 
the Bianchi I model are encoded in three triad components pi and momentum components c$, related to the metric 
components as 

Pi = a 2 a 3 , p 2 = aia 3 , p 3 = a\a 2 , Ci = jai , (17) 

where 7 is the real- valued BarberoTmmirzi parameter and represents an ambiguity parameter of loop quantum gravity. 
Black hole entropy calculations can be used to fix its value. In terms of these variables, the Hamiltonian is given by 

H = - — 9 \_ (c 2 P2C3P3+ClPlC 3 p 3 +C 1 p 1 C 2 P2) + Hm, (18) 

8nGYy/pip 2 p3 

where TIm is the matter contribution to the Hamiltonian. Einstein's equations can then be derived from Hamilton's 
equations, which explicitly for this system are 

Pi = SttG^ , b % = SirGj^ . (19) 
oci dpi 

The Hamiltonian must also vanish for the system: 

n = o. (20) 

The Hamiltonian for the matter contribution is proportional to the energy density of the matter, so the magnetic 
Hamiltonian is given by 

n B = a 3 p B = ^. (21) 

We will also consider a perfect fluid with constant equation of state w. This can be added to the matter Hamiltonian 
by first solving the conservation equation to give 

p = Ca-^ 1+W ^ , Hfiuid - Co- 3 ™ , (22) 

where C is a constant. 

Analyzing the system at the level of the quantum difference equations of LQC for this model would be highly 
challenging, given the complexity of the Bianchi I equations. We thus consider approximate semi-classical equations 
of motion that incorporate loop quantum modifications. These effective equations have been shown to be very good 
approximations for the case of isotropic cosmologies sourced by a massless scalar field 0, H, 0, H, H , and the results 
have been extrapolated to more complicated models. The corrections modify the general relativistic Hamiltonian (118|) 
to be of the form [l(| 

1 rsin(/22C2)sin(/2 3 C3) 1 
Heft = - „ 5 _ < — p 2 P3 + cyclic terms } + H M - (23) 

The parameters \ii are the key ingredients determining the quantum corrections. It is easy to see that in the limit 
[li — > 0, the classical Hamiltonian (|18p is recovered. The pi parameters are assumed to be functions of the triad 
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components vj. and their precise specification is an ambiguity of the quantization. Two possible constructions are 
discussed in [lCj , although one of them has been argued to be more physical on the grounds of certain gauge invariance 
considerations in [ll| and of a more rigorous construction of the quantum theory in [l2j ■ We will focus on that scheme 
in this paper. The particular form is 



Mi 



(24) 



where A is a constant that is typically related to the minimum area gap of loop quantum gravity. In this paper we 
assume that A = O(l) in Planck units; the precise value will not affect the qualitative results. The effective equations 
of motion can be derived as in the general relativistic case, using Hamilton's equations (|19p and the vanishing of the 
Hamiltonian. 




FIG. 1: The pure-magnetic axisymmetric case. Top: The scale factors are shown on the left, and the volume factor V — axa^ao, 
is shown on the right. Dashed lines indicate the classical behavior, a\ — > 0, a-2 = az — > const, and V —* 0. Solid lines show the 
effective loop quantum solutions. Quantum effects regulate the singularity leading to a bounce in 02, 13 and overall expansion 
of the universe. Note that a 1 continues to decrease after the bounce. Bottom: The energy density as a fraction of critical (left) 
and the shear energy density E 2 = a^ v a^ v /2 (right). This shows the non-conservation of shear anisotropy through the bounce. 

The equations of motion are sufficiently complicated to not allow for an analytic solution. Despite that, some general 
conclusions can be made from the form of the equations. First, as shown in [id ], the vanishing of the Hamiltonian (|23|) 
immediately implies a bound on the energy density of the matter. This arises from the bound in the sin terms of the 
constraint. The precise bound is the same critical density that characterizes the bounce in the isotropic models: 

P ^sJfA- (25) 

The total energy density of the matter (magnetic plus fluid) must be below this value. This is an indication that the 
classical singularity (where the energy density diverges) is removed and replaced by a bounce. The second conclusion 
from the effective equations, is that if the matter has zero anisotropic stress, the shear term is conserved before and 
after the bounce. Since the magnetic field has non-zero anisotropic stress, Eq. ([9]), this behavior is not guaranteed. 
In the next section, we present numerical results. 



IV. NUMERICAL SOLUTION OF THE EFFECTIVE LOOP QUANTUM EQUATIONS 



The first case is an axisymmetric spacetime (02 = 03) sourced only by the magnetic field. Classically, the singularity 
is characterized by Eq. (|15p . We use initial conditions corresponding to a classically collapsing universe approaching 
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FIG. 2: The pure-magnetic non-axisymmetric case. Dashed lines show the classical behavior, a±,a2 — > 0, 03 — » 00 and — > 0, 
and solid lines show the effective loop quantum solutions. 

the classical singularity. The solution is shown in Fig. [1] The quantum solution matches the classical well until 
the singularity is approached. Then the quantum effects act repulsively preventing a± from reaching zero, and 
leading to bounces in 02 = 03. The volume factor confirms that a bounce replaces the classical singularity. The 
post-bounce expanding universe has 02 = 03 — * 00 and a\ — > 0, as in the classical case. Thus the quantum effects 
join a classical contracting branch with an expanding classical branch. The energy density shown in Fig. [1] remains 
bounded below the classical critical density p c as expected from analytical considerations. Finally, the shear energy 
density S 2 = er^er^ /2 is shown to remain finite through the evolution, but is not conserved through the bounce. This 
is in contrast to the pure-fluid case [Iol |. where shear is conserved. The difference arises from the non-zero magnetic 
anisotropic stress, Eq. ([5]), which leads to production of shear anisotropy. 
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FIG. 3: As in Fig. [2] but with different, Kasner-like initial conditions. 

Figure [2] shows the solution for the pure-magnetic non-axisymmetric case. The classical singularity is described in 
Eq. (fl6|) . Once again a bounce occurs in the volume near the point of the classical singularity. With the choice of 
initial conditions in Fig. [21 classically a 3 — > 00 while 02 — > at the singularity. With the quantum effects, a\ and a-i 
are repelled from zero, and a-i bounces. As in the axisymmetric case, the post-bounce regime is an expanding universe 
with ci2, 0,3 expanding and a\ contracting. The overall behavior is qualitatively similar to the axisymmetric case, and 
again the energy density of the magnetic field is bounded below p c and the shear is not conserved. 

An alternative non-axisymmetric choice of initial conditions is shown in Fig. [3] The initial conditions are chosen 
to be in a Kasner phase, where two directions are contracting and one is expanding, while overall there is contraction 
in the volume. The qualitative behavior is qualitatively similar to the first non-axisymmetric example in Fig. [21 

As a final case, we included a dust perfect fluid, w = 0, with the magnetic field. We chose non-axisymmetric initial 
conditions with the volume collapsing. Figure [3] shows a bounce qualitatively similar to the pure-magnetic case. The 
shear term remains finite, but is not conserved. At late times in the post-bounce expansion phase, the dust begins to 
dominate the evolution and the universe isotropizes, since the ratio of the separate Hubble rates tends to one at late 
times. 
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FIG. 4: Magnetic field and dust, non-axisymmetric. Top: The classical singularity is avoided via a quantum bounce (left), 
and the shear is not conserved through the bounce (right). Middle and bottom: The ratios of expansion rates, showing the 
late-time, post-bounce isotropization due to the dust. 

V. CONCLUSIONS 

We have extended the effective LQC treatment of Bianchi I cosmologies by including a homogeneous magnetic 
field. We have studied dynamics using the approximate effective equations of motion that capture features of the 
(as yet unknown) true quantum LQC dynamics. Thus our results are approximate, and a more rigorous quantum 
construction would be needed to fully validate them or provide additional corrections. 

The effective equations indicate that the singularity-avoiding bounce is not spoiled by the inclusion of a homogeneous 
magnetic field. Extending the results of [T(|], we showed that shear anisotropy does not blow up as the classical 
singularity is approached, but remains finite through the entire evolution. In contrast to the pure-fluid case [Iol |. we 
showed shear is no longer conserved through the bounce, due to the anisotropic stress carried by the magnetic field. 
Our results indicate an interesting evolution of shear, possibly with a net generation of shear, but further study is 
needed to check whether this is generic. When we add a dust fluid to the magnetic field, the qualitative behavior 
through the bounce is unchanged. However, at late times after the bounce, the dominance of the dust ensures that 
the universe isotropizes, unlike the pure-magnetic case. 
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